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ABSTRACT 


Two  problems  relevant  to  the  magnetotelluric  detection 
of  high  altitude  nuclear  detonations  are  treated  in  this  report: 

(1)  A  computer  program  is  described  for  solving  the  poloidal 
magnetohydrodynamlc  wave  equation  with  a  source  current  present 

in  a  realistic  though  axisymmetric  model  of  the  outer  ionosphere. 
Arguments  are  given  that  the  poloidal  mode  of  the  axisymmetric 
model  should  produce  the  relevant  features  of  a  more  realistic 
asymmetric  model. 

(2)  The  effects  of  collisions  on  the  propagation  of  the 
wave  signal  through  the  lower  ionosphere  is  studied  experimentally 
in  a  laboratory  plasma.  The  facilities  constructed  for  the 
purpose  and  the  experimental  program  are  described  in  this  report. 
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A.  Theoretical  Study  of  Magnetohvdrodynamlc  Wave  Propagation 


1.  Introduction 


This  work  has  been  sponsored  by  the  Advanced  Research 
Projects  Agency  under  the  Nuclear  Test  Detection  Program.  The 
magnetotelluric  detection  of  high  altitude  nuclear  detonations  is 
far  less  advanced  than  many  other  detection  methods.  One  of  the 
primary  reasons  for  this  is  that  a  complete  theoretical  understanding 
of  the  propagation  of  a  magnetic  disturbance  in  the  ionosphere  is 
lacking. 

This  is  the  subject  of  the  present  theoretical  study. 

We  are  concerned  with  the  radiation  of  magnetohydrodynamlc  waves  by 
a  current  source  located  within  the  ionospheric  cavity  surrounding 
the  earth. 

In  particular,  we  treat  only  the  poloidal  or  isotropic 
mode  which  is  uncoupled  from  the  toroidal  or  guided  mode  for  an 
axlsymmetric  geometry.  This  simplification  originates,  of  course, 
in  the  Interests  of  tractability,  but  some  recent  studies,  described 
in  section  2,  indicate  that  the  amplitude  of  the  guided  mode  is 
generally  several  orders  of  magnitude  smaller  than  that  of  the 
isotropic  mode  if  the  ionospheric  medium  is  not  too  highly  asymmetric. 

The  problem  we  are  treating, then,  is  an  ordinary  waveguide 
problem,  but  for  a  realistic  variation  of  the  ionospheric  parameters 
(ambient  ion  density  and  magnetic  field  strength)  it  is  still  not 
amenable  to  analytic  solution  and  the  wave  equation  must  be  solved 
numerically.  This  computation  using  the  method  of  lines  (straight- 
line  method),  attributed  to  M.  G.  Slobodyanskii^  ,  is  described  in 
section  4. 
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First,  however,  in  section  3  we  solve  the  poloidal  wave 
equation  analytically  for  the  special  (and  unrealistic)  case  in 
which  the  Alfven  speed  varies  linearly  with  th»  distance  from  the 
center  of  the  earth.  The  object  here  is  twofold:  (1)  We  can  get 
some  assurance  that  the  computer  program  described  in  section  4  is 
working  correctly  by  testing  it  on  this  special  case.  (2)  We  can 
get  some  feeling  for  the  way  in  which  the  orders  of  magnitude  of 
the  various  ionospheric  parameters  affect  the  wave  propagation. 

In  particular,  the  analytic  works  of  section  3  takes  the  finite 
ionospheric  conductivity  into  account. 

Finally,  in  section  5  the  analytic  and  computer  studies 
are  compared  for  the  special  case  previously  mentioned  and  the 
computer  results  for  more  realistic  ionospheric  models  are  described 
and  Interpreted. 
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2^__jome_Recent_DeveLopments  In  the  Theory  of  Magnetohydrodynamic 
Wave  Propagation  . ----  ---  - 


The  calculational  work  has  dealt  with  solving  the  axisym- 
metric  pololdal  equation.  While  consideration  of  this  special  case 
certainly  does  not  yield  a  complete  description  of  magnetohydro¬ 
dynamic  wave  propagation,  there  are  good  arguments  as  to  why  this 
case  is  an  important  element  of  any  such  description. 

The  basic  magnetohydrodynamic  equations  are  given  in 
section  3.  If  the  assumption  is  made  that  the  field  amplitudes 
are  independent  of  uf  ,  then  the  magnetohydrodynamic  wave  equation 
breaks  up  into  two  independent  equations.  One  of  these  describes 
the  pololdal  mode  which  will  be  treated  in  more  detail  below.  The 
second  is  the  toroidal  mode  which  has  properties  making  it  similar 
to  the  guided  mode  for  plane  wave  propagation  in  a  constant  density 
medium  with  infinite  conductivity  and  a  constant  magnetic  field. 

For  the  toroidal  mode  in  a  dipole  field,  the  equations  are 


To  -*■  jjk  •  V  -  Ly 


Hirpn  dtfj 
dt 


and 


-  F  AT 

vt  L  -  J  ? 


(2.1) 


(2.2) 


Since  £»V  is  the  only  operator  present  in  these  equations, 

0  2 
a  change  of  variable  is  suggested  from  r ,0  to  v,0  where  v  =  sin  0/r. 

With  (i  —  cos  0,  these  equations  become 


^  ,  (H  ^ **“  ^ 


(2.3) 


and 


**V(  f 


H 


#1 


T.) 


W  Am4*/. 


(2.4) 


Defining  g(ji,v)  «  v^/sin^e  and  02(n,v)  -  4*p  ( 1-ji2)  *7(Mv^) 2  the 

f  o 

toroidal  equation  has  the  form 


-  ft  T$ 

3t 


(2.5) 


From  this  equation,  it  is  easily  seen  that  the  eigenvalues  arising 
in  the  solution  of  any  boundary  value  problem  are  functions  of  v, 
denotes  the  second  derivative  of  g  with  respect  to 


since 


|i  with  v  being  held  constant.  Since  v  s  sin  0/r  is  a  function  of  the 
spatial  coordinates  and  is  essentially  a  parameter  which  labels 
different  geomagnetic  lines  of  force,  we  have  the  unusual  situation 
that,  for  toroidal  axisymmetric  wave  disturbances,  the  eigenvalues  are 
a  function  of  position. 
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(o) 

This  condition  has  been  discussed  by  Radcski  and  McClay  '  . 

In  carrying  out  an  analysis  using  a  simplified  geometry,  they  note 
that  the  charge  density  and  the  current  density  parallel  to  the 
magnetic  field  lines  are  both  proportional  to  derivatives  with 
respect  to  v  of  field  quantities.  Now  the  field  quantities  vary  in 
time  as  e*^  where  ■  tu(v)  from  Eq.  (2.5).  Thus  we  have  a  situation 
where  the  charge  and  parallel  current  densities  vary  as  5®  ,  that  is, 

they  grow  indefinitely  with  time.  This  contradicts,  among  other 

icot 

things,  the  initial  assumption  that  all  quantities  vary  as  e 
Thus  one  must  conclude  that  under  the  hydromagnetic  approximation, 
the  toroidal  field  must  be  zero.  Stated  differently,  in  order  for 
the  eigenvalues  of  the  toroidal  mode  to  depend  on  the  cavity  size, 
rather  than  on  the  label  v  of  a  particular  magnetic  field  line,  it 
is  necessary  to  introduce  into  Eq.  (2.5)  derivatives  with  respect 
to  v,  so  that  the  toroidal  field  amplitude  is  coupled  across  field 
lines.  Once  this  is  done,  however,  the  simple  model  of  Alfven  waves 
propagating  along  geomagnetic  field  lines  and  perhaps  bouncing  back 
and  forth  between  conjugate  points  with  a  characteristic  bounce 
period  must  also  be  laid  to  rest. 

One  means  of  introducing  derivatives  with  respect  to  v  in 
(2.5)  is  to  allow  the  presence  of  a  finite  electric  field  parallel 
to  the  magnetic  lines  of  force,  that  is,  change  the  simple  equation 

L  v  x  & c 

c  ^  ^  (2.6) 


to  allow  for  other  effects  such  86  Hall  effect,  electron  pressure, 

finite  conductivity,  etc.  If  the  Hall  effect  is  introduced,  so  that 

(3) 

a  parallel  electric  field  component  is  introduced,  it  may  be  shown 
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that  the  poloidal  field  is  not  changed,  but  that  the  toroidal  fie..d 
is  coupled  to  the  poloidal  field  and  that  the  toroidal  field  vanishes 
in  the  limit  of  the  hydromagnetic  approximation.  Also  for  geophysical 
parameters,  the  toroidal  field  is  small  compared  with  the  poloidal 
field. 

Thus  we  are  led  to  conclude  for  axisymmetric  magnetohydro¬ 
dynamic  wave  propagation  that  it  is  the  poloidal  mode  which  is 
important. 

The  axisymmetric  poloidal  and  toroidal  modes  are  the  m  *  0 
modes  which  result  when  the  field  amplitudes  are  assumed  to  vary  as 
eim^  (where  m  is  an  integer  and  cf  is  the  third  coordinate  of  the 
usual  (r ,6 spherical  coordinate  system).  The  higher  order  modes 
such  as  m  *  1,2, etc.  cannot  be  so  easily  neglected.  However,  one 
can  make  two  statements  about  the  higher  order  modes: 

1.  Generally  speaking,  unless  a  particular  symmetry  is 
specifically  picked  out  by  the  nature  of  the  physical  problem, 
modes  with  smaller  values  of  m  will  be  more  important  than  modes  with 
larger  values  of  m. 

2.  Starting  with  the  magnet  ohydrodyna.nic  wave  equation 
for  propagation  in  a  dipole  field,  Dungey  has  shown  that  in  the 
limit  of  large  values  of  m  (i.e.,  for  high  azimuthal  asymmetry) 
there  is  only  one  mode  propagating.  This  mode  has  the  same  field 
variables  describing  it  as  does  the  m  -  0  poloidal  mode.  On  the 
other  hand,  as  already  indicated,  for  m  large  the  field  amplitudes 
which  described  the  toroidal  mode  for  m  *  0  are  small.  A  further 
very  important  property  of  the  surviving  mode  in  the  large  m  case  is 
that  its  wave  equation  has  a  form  quite  similar  to  Eq.  (2.5).  In 
particular,  the  wave  amplitude  in  the  highly  asymmetric  case  do  not 
depend,  to  a  first  approximation,  on  derivatives  perpendicular  to 
the  magnetic  field  lines.  Thus,  the  eigenfrequencies  depend  upon 


-6- 


spatial  coordinates  as  before.  Since  this  leads  to  an  unstable 

growth  of  the  mode  we  are  led  to  conclude  that  these  modes  have 

(3) 

near-zero  amplitude  in  a  first  approximation  . 

From  this  we  are  led  to  believe  that  the  m  =  0  poloidal 
mode  whose  equation  we  have  solved  numerically,  and  the  other  low 
order  modes  are  dominant  in  problems  of  magnetohydrodynamic  waves 
propagation  in  a  dipole  field  in  he  vicinity  of  a  spherical  earth. 
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3.  Analytic  Solution  of  the  Poloidal  Magnetohvdrodvnamic  Wave  Equation 
with  Finite  Atmospheric  Conductivity  for  a  Special  Case 


The  basic  equations  are  as  given  in  a  previous  report 


(6) 


with  the  addition  of  a  finite  scalar  conductivity. 


(3-1) 


(3-2) 


(3-3) 


(3-4) 


From  (3-1)  and  (3-2)  : 


or  using  (3-4)  and  (3-3) 


B  •  X  Ae  /  vx  Vxt  -b  8,  /  %0  Xwr 

A-  **  A*  - 


jJ 


--t/rf.  ii./ts'le)  =‘ 


(3-6) 
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but  from  (3-5) 


n  =--*2  Ft?*? 

~7tr  7*  f-r 


(3-7) 


so  that  (3-6)  becomes: 

&oX  &o  xyx7*£  't'leXAc  *2 

A  A.  A»  <V  .V  .  £  - 

O B 

-  17  f,  **£  (3-8) 

0~  ji 

-ir‘f  °l~  +  ^  v*y*<>i  ) 

'  ft-e  a  fr  jitr  j 


loot 

Taking  an  e  time  dependence  this  becomes: 


£&>  *  ~ '%  c  ?*v*€ 


(3-9) 


4-  4-77 f0  *0  t/F  =  -[  tlo  */&  cicoJ  trf£V& 

cr  c.L 


We  take  J  to  have  only  a  component  and  B  to  have  *r  and 

~  r  ~o 

components.  All  quantities  are  assumed  ^  independent. 
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where: 

/3  -  ff if  £i  7 

"If*  (3-i3) 

-2 

As  a  special  case  we  assume  V^e(  r  and  oe( r  so  that  p  is 
a  constant* When  J  ■  0  the  homogeneous  equation  is  separable  and  the 
orthonormal  radial  and  angular  eigenfunctions  satisfying  the  boundary 
conditions  E  ■  0;  r  ■  r^rji  ®  are  found  to  be: 


-10- 


14,-H 


%  v  iqtH) 
hffCy)-  )  l/i-J 


tic*) 


^  r^i 


'vMvihrMX  1 


The  corresponding  eigenvalues  are: 


We  take  a  current  density: 

J  =  -ZT  gTCt~*>)  J~Cx-x9) 
Br*0i 


where  I  is  the  total  current. 

The  solution  to  (3-12)  is  then  just  the  Green's 


-Zi?zf  te(*)4«crj 

c  6  **S/  Atm 


(3-14) 


(3-15) 


(3-16) 


(3-17) 


function: 


(3-18) 
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t*''  t €/4#J  j£+^J)Z+#*s>~/*] 


t^R_^  J  AfpJUjr  ^ 

j£  '  U«j  '  J 


The  l  summation  can  be  done  with  the  aid  of  the  following 
formula  from  Magnus  and  Oberhetinger^  (the  m  sum  could  also  be  done 
analytically  but  this  does  not  prove  to  be  useful): 


/%'  Or)  fCw  ■<> 

*  f  (-0  *f'  -(?*■**//' J 1%  '(Xjt  'r>0 


(3-19) 


In  equation  (3-18)  we  put  =  f- 

so  the  sum  to  be  done  is: 

y  2. 4*/  /*£  CAe)^Of) 

-  2  (s-t+i) (-ij  *  'Ob)  ft'c-x) 

*='  W&tv 

where  we  have  used  the  relations 

/¥  <x)  ^-(-0  </%'(-*) 

/%'(-*)  -  eft*)  &'c~*) 


-/>*#:  i 


(3-20) 


(3-21) 
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(3-22) 


In  equation  (3-19)  we  put  v  ■  = 

fr  )  P~'  (-X )^-  T~/'±Ut  Tt  * 

i  hj  (  fc&'J  ~'c-x J 

*=/ 

(3-23) 


comparing 


f  (s-e+d  ft'(Xo)  Z'OrJ 

if 


^  T 


“»)  C,  \ .  (n) 


w 

'%*'•**' 


-  -V _  P  /  (y  )  P' 


C-x) 


~(«m‘  -£_) 


where  we  have  again  used  (3-22). 

Thus  we  finally  obtain  Crom  (3-18): 


(3-24) 


zgz(A  */(*y#z 

M  *  /  ^ 


4*+v/  b#’  lu  jS  I  At-  >2. 

I  *  /  / 


(3-25) 


—  v> 

vAii 


; 
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The  i  summation  in  (3-18)  could  have  been  done  directly  by 
writing  the  sum  as  a  contour  integral  and  using  the  Watson  transforma¬ 
tion.^  The  sum  in  (3-25)  would  then  be  the  so  called  residue  series 

We  write  p  ■  pr  +  ip1  and 


V  1  ^ 

-  (hJ+f.  V  V  Z  t  -  74 

-  '//*£■ 

2. 


(3-27) 


7^ 


Similarly  we  find: 

«C  =  ^  +(■]**. 

(3-28) 
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is  large. 


Now  /  /  is  generally  a  large  number  because  X  •  — 

This  might  not  be  true  for  one  value  of  m; 


ro) 


(  y-  r/S  +77~<r  Igl 

this  is  a  special  case  which  can  be  treated  separately. 
For  large: 

4*  u°)  * 

/^W/ J 


^  ol 


m 


i/z  i/mr  ) 


,  but 


(3-29) 


u>  c 

The  situation  of  interest  is  —j-  «  1  and  since  we  are 

\ 

excluding  m  too  close  to  In  ^2.  we  have  a  «  b  so  that  (3-27)  and 

rl  X 


(3-28)  become: 


.  >  ^  JL.  £5  3  T  *-Mur  tz7\ 

*} 


«m  ~^z  *  yi[/  -  \  1 J  '/z 


)**  *- 


*i 


(3-30) 


(3-31) 


>? 


$ 


and: 
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(3-32) 


If  (  )  >  1  then  (3-29)  gives: 

l ffe-  e 


1  _? -  -  g-**"  {€>-<%) 

^~<*nt  ILuzq  •*&& 

r  ^  s<n,] 

*  c^«t»T c r-e)  -/  -*<S  (r-g)] 

'AjSSn&s 


(3-33) 


(3-34) 


(3-35) 


&1*)  +'*#**<* 
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Clearly,  the  only  terms  of  importance  for  6  not  too  close 

* 

to  in  (3-25)  are  those  with  the  smallest  .  From  (3-30) 

and  (3-32)  this  implies  /  K 


fe)  <1‘ 

'  r,  ' 


Thus  the  summation  in  (3-25)  may  be  truncated  at  an  integer  m  where  m 
X  r 

is  the  order  of  —  In  2  . 

n 


Putting  (3-35)  in  (3-25)  we  obtain: 


.  •  .  r 

assuming  1  «  %  « 

The  magnetic  perturbation  at  the  lover  cavity  boundary  is: 


(3-37) 


ItSeFl'lS  I  --  >J:  L  *  (yeJ! 

'  -  M J  ~  lr-_Yt  To  Y  <>r[  'r-T 

x'1' i.  r,~** 


ftl*/ 


4  \r 


>Hrt  twrj?'  f  ) 
V  Hit  '' 


1 
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For  a  low  altitude  eource  fry  / 


K 


&&  ~  ^  -2^-  li y  >;  ^t-'r  )  (^<+^  +4m,&c  j  ''f'z 

/  *  z  <L5>^  £  ^  j  *  ;** u£«-«)J 

**"1^  '  (3-38) 

-  $L  jz  (tn  **■ )  '  KT^  V  (*<>  j(u*LtT44#i^j-'/l 


mt  , 

P  / 


-/  (T7ij^ 

^  ICs-^j'^-^j  s  X 


The  impediment  to  summing  this  series  is  the  cos  factor  for 

2  — 
while  1  -  x  may  not  change  appreciably  as  m  goes  from  1  to  m, 

2  h 

X(1  -  x  )  will.  Finally,  we  note  that  if  cr  is  small  enough  the  m  -  1 
term  dominates  through  the  exponential. 

We  defer  further  consideration  of  the  analytic  solution  until 
after  the  computer  solutions  have  been  treated  in  the  next  section.  We 
will  then  c?poare  the  analytic  and  computer  solutions  and  indicate 
how  the  orders  of  magnitude  of  the  actual  ionospheric  variables  are 
expected  to  affect  magnctohydrodynamic  wave  propagation. 
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For  a  low  altitude  source 


H7 


) 


w 

^  Z’Q*.  Vjjy  K  ('6-y  J ***■%, 

/  *  z  ^  <<&-* ;  f  ;** 

/r7  (3-38) 

-  2ljz  (tn  *£  )  1  r*''*.  £  ”$£  (fe  -)  /J£u*ut44*qf  j 

,  e  ’  #»  *  ^  V  V* 


C~ 

fir 

/>  j 

n) 


-.,  C*iJ 


)  r/  ^  -t% j  ."r  = 


=./JwZZ 

At  >i 


The  impediment  to  summing  this  series  is  the  cos  factor  for 

2  - 
while  1  -  x  may  not  change  appreciably  as  m  goes  from  1  to  m, 

2  k 

X(1  -  x  )  “  will.  Finally,  we  note  that  if  cr  is  small  enough  the  m  ■  1 
term  dominates  through  the  exponential. 

He  defer  further  consideration  of  the  analytic  solution  until 
after  the  computer  solutions  have  been  treated  in  the  next  section.  We 
will  then  compare  the  analytic  and  computer  solutions  and  indicate 
how  the  orders  of  magnitude  of  the  actual  ionospheric  variables  are 
expected  to  affect  magnetohydrodynamic  wave  propagation. 
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4.  Numerical  Solution  of  the  Wave  Equation 


The  computer  program  which  we  have  developed  is  directed 
towards  solving  a  simple  model  incorporating  several  properties 
of  the  actual  ionospheric  environment:  (1)  Spherical  geometry 
(2)  Dipole  Magnetic  Field. 

For  a  realistic  treatment  of  magnetohydrodynamic  wave 
propagation,  it  is  not  sufficient  to  use  ray  tracing  procedures  or 
a  plane-layered-medium  approximation.  This  is  easily  seen  since 
the  wavelength  A  ■  V  /f  where  V  is  the  Alfven  velocity  and  f  is 

fl  8 

the  frequency  is  the  order  of  1000  km  or  greater  for 
8  . 

V  »  10  cm/sec  and  f  «  1  cps  or  less.  This  is  the  scale  over 

3 

which  the  geomagnetic  field  varies  as  well  as  the  distance  for 
which  earth-curvature  effects  become  important.  Thus,  not  only  is 
it  important,  if  one  wishes  to  be  realistic,  to  consider  the 
earth's  curvature  and  the  dipole  nature  of  the  geomagnetic  field, 
but  the  magnetohydrodynamic  equation  must  be  solved  without  using 
a  geometrical-optics  approximation. 

The  basic  equation  is  (3-8)  with  o  *qp: 


KAkHW'Z)]}  =  —  [g  *(>0> 

W  ^  4  >  /0ca  ^  zt'l 


(4-1) 


The  geomagnetic  field  is  taken  to  be  of  the  form 


The  equation  separates  into  two  independent  equations  if  d  and  £  are 

assumed  to  be  independent  of  (i.e.,  assuming  axial  symmetry). 

One  of  these  equations,  for  the  toroidal  mode,  has  already  been 

2  2 

discussed.  The  other  equation  is  (3-12)  with  (J  =  K  (r)  r  where 
K2(r)  =  4 «!Pq  (r)  r V  : 

M2 


(4-3) 


This  equation  has  been  solved  numerically  using  the  method  of 
lines.  In  this  method  the  equation  is  gridded  along  linos  of 
constant  x.  Thus,  with  E^r.x)  *  E^(r),  we  have  a  set  of  equations 
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where  1  <  j  <  N  +  1  and  N  ■  2/h.  h  is  the  step  size  in  x  in  the 
angular  direction.  The  boundary  conditions  E^(r)  ■  E^^(r)  ■  0 
are  imposed  since  E  ^  is  zero  at  9  ■  0°  and  180°.  The  set  of 
equations  (4-4)  are  coupled  second  order  differential  equations  which 
must  be  integrated  in  r.  The  boundary  conditions  E^(r^)  ■  E^r^  ■  0 


are  imposed  at  the  inner  and  outer  boundaries.  The  source  current  J 
we  have  been  using  is  zero  everywhere  except  on  a  line  of  constant  x. 
C.i  this  line  it  is  zero  except  between  the  radii  rfl  and  r^,  where  it 
is  constant.  Other  sources  may,  of  course,  be  used. 


j 
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The  basic  problem  of  solving  a  set  of  coupled,  ordinary 
differential  equations  for  a  two  point  problem  can  be  attacked  in 
several  different  ways.  The  method  used  in  this  report  is  discussed 
in  the  Appendix. 

Many  of  the  cases  run  up  to  the  present  involve  taking 

2 

Kr  ■  constant  ■  X  and  ignoring  the  l+3x  (i.e.,  taking  it  equal  to  1) 
in  Eq.  (4-4) .  The  reason  for  this  particular  choice,  as  already 
pointed  out,  is  that  the  equation  is  separable  and  can  be  solved 
analytically.  A  comparison  between  the  analytical  and  the  numerical 
solution  is  made  in  the  next  section. 

Also,  after  the  computer  program,  as  described,  was 

running  an  entirely  different  program  was  set  up,  similar  to 

(8)  f9) 

Greenfield's  ■  and  to  Chapman  and  Bostick's k  in  which  the  partial 

differential  equation  was  differenced  in  two  dimensions.  A  compa’-i- 

son  of  the  results  of  this  program  with  the  one  using  the  method  of 

lines  gave  excellent  agreement.  This  is  another  check  on  the 

correctness  of  the  program. 
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5.  Numerical  Results  and  Interpretations 


We  first  consider  the  case  in  which  X  =  ru>  is  a  constant  and 

VA 

the  conductivity  is  assumed  to  be  infinite.  The  bounding  radii  of  the 
cavity  are  6,520  km  (150  km  altitude  taking  rg  *  6370  km)  and  19,880  km. 

Throughout  this  section  the  source  current  is  located  on  the  line 
x  *  cos  0  *  -.6  (126°  53')  between  6,674  and  6,857  km  (304  to  487  km 
altitude) . 

Figures  (5-1)  through  (5-4)  show  the  results  for  \  ■  10  and 
(5-5)  through  (5-8)  for  X  =  100.  The  scale  of  the  amplitude  is 
arbitrary  and  the  field  is  plotted  as  a  function  of  angle  for  two 
values  of  constant  radius  and  as  a  function  of  radius  for  two  values  of 
constant  angle.  The  remarks  which  follow  pertain  as  well  to  the  field 
amplitudes  for  other  values  of  radius  and  angle  which  have  not  been 
plotted  although  they  are  contained  in  the  computer  output. 

Unfortunately,  obtaining  the  field  as  a  function  of  angle 
is  relatively  time  consuming  since  the  time  varies  as  the  square  of 
the  number  of  angular  points  and  only  linearly  with  the  number  of 
radial  points.  This  is  the  resul  of  gridding  along  lines  of  constant  x 
rather  than  r  as  described  in  section  4.  The  following  figures  there¬ 
fore  generally  give  only  a  sketchy  indication  of  the  angular  variation 
although  the  radial  variation  is  given  in  some  detail. 

Since  the  wavelength  is  smaller  for  X  *  100  than  for  X  ■  10 
it  i  'i  not  surprising  that  the  X  ■  100  curves  are  more  oscillatory. 

The  most  striking  feature  of  these  plots  is  that  the  field  is  nearly 
separable.  For  example  the  number  of  nodes  in  the  radial  function 
is  independent  of  the  value  of  the  angle.  Separability  generally 
seems  to  be  poorest  in  the  vicinity  of  the  source. 
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r  -  8,153  km 
5  angular  points 


J 


r  =  16,963  km 
5  angular  points 


This  observation  suggests  that  one  term  in  the  series  (3-25) 
is  dominant,  and  this  is  in  fact  the  case.  The  dominant  term  is, 
roughly  speaking,  the  one  with  the  smallest  value  of  !°*/»/  and  its 
dominance  apparently  arises  through  the  factor  occuring 

in  each  term  of  (3-25).  The  other  terms  seem  to  be  of  secondary 
importance  because  they  are  all  of  the  same  size  but  of  random  phase. 


Obviously,  this  can  only  be  true  if  we  are  not  near  an 
eigenfrequency,  for  then  o fa  where  &  is  an  integer,  and 

the  dominance  of  one  term  occurs  through  the  (c*0<-  0^ /pj  factor 
in  each  term  of  (3-25) . 


The  actual  dominant  m  value  (1R)  may  be  such  that  cx.  is 
either  the  smallest  real  or  the  smallest  imaginary  value  of  .  We 
have  at  present  no  way  of  predicting  which  of  these  two  possibilities 
will  occur. 


Figures  (5-9)  through  (5-12)  show  two  cases  where  —  is  the 
smallest  imaginary  value  of  . 


In  the  first  case  X  -  90  and  is  again  19,880  km  in  the 
second  case  X  ■  59.5  and  ^  is  decreased  to  8,370  km.  The  major 
difference  between  these  figures  and  Figures  (5-1)  through  (5-8) 
is  that  the  angular  function  is  now  damped  going  away  from  the  source. 
This  makes  sense  because  when  o(m  is  imaginary  the  angular  variation 
of  the  mth  term  of  (3-25)  is  the  conical  function 

which  is  monotonic  decreasing.  A  few  of  these  functions  are  tabulated 
in  our  previous  report^ for  small  /°^n// 

/  1 

For  /<*„/  large  the  factor 

which  gives  the  angular  variation  of  the  mth  term  is  just  proportional 
to  (&  "Go) 


The  fall  off  in  Figure  (5-12) 


VXcZe 
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mm 


should  in  theory  be  characterized  by  2.C  »  fact>  it 

is  more  like  /**}/  “i  /O  .  This  discrepency,  we  think,  arises 
from  the  difficulty  the  difference  equations  described  in  section  4 
have  in  following  the  rapidly  decreasing  exponential  function. 

Table  (5-13)  illustrates  the  results  of  a  number  of  runs 
with  r2  *  19,880  km  and  with  various  \  values. 

The  number  of  npdes  in  the  radial  function  (second  column) 
is  determined  from  the  computer  cutput.  is  found  from 

(3-31)  and  (3-32)  with  mi  equal  tne  number  of  nodes  listed  in  the 
preceding  column  plus  one.  In  every  case  with  constant  \  investigated 
thus  far  has  been  either  the  smallest  real  value  of 

or  the  smallest  imaginary  value. 


In  the  last  column,  the  number  of  nodes  anticipated  in 
the  angular  function  is  computed  by  eisuming  that  when  is 

real  the  angular  function  will  behave  like  f C~x)  where  i  is 

the  nearest  integer  to  —  V-  nt 

In  only  one  case  have  we  used  enough  angular  points  to 
adequately  test  this  hypothesis.  For  this  case  in  which  20  lines  of 
constant  x  were  used  and  the  parameters  were:  \  ■  59.5  r2  *  19,880, 

we  find  °%  -6*5  and  i  ■  6.  The  angular  function  does,  in  fact, 

have  5  internal  nodes,  in  agreement  with  our  hypothesis. 


If  it  is  true  that  the  mth  term  alone  of  (3-25)  is  a 
reasonably  good  approximation  to  the  field  we  have: 

r/j 


E  * 


l /'(*  '■ 

'  Co*, 


(5-1) 
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TABLE  (5-13) 


19,880  km 


r 


2 


Number  of  Nodes  of 
the  Radial  Function 
Excluding  the  End 
Points 

—  is  the 

smallest  - — 
value  of 

xm 

i  Number  of  Nodes 
of  the  Angular 
Function  Exclud¬ 
ing  the  End  Points 

2 

real 

4 

6 

real 

10 

10 

imaginary 

damped 

13 

real 

5 

16 

real 

16 

20 

real 

5 

20 

real 

9 

24 

imaginary 

damped 

26 

real 

23 

31 

imaginary 

damped 

34 

real 

15 

The  addition  of  finite  conductivity  according  to  (3-30) 
just  adds  a  damping  factor: 

CO  — 

e~  I**) 

(5-2) 


to  (5-1)  if  o  is  not  too  small  and  if  ’>  > 

We  turn  now  to  the  case  which  introduces  the  angular 

variation  characteristic  of  the  dipole  field:  X  *  Xq  , 

'  2 
1+JX 

where  X  is  constant, 
o 

’’’he  results  are  shown  in  Figures  (5-14)  through  (5-20) 
for  Xq  *  59. j  and  r^  ■  19,880,  The  primary  effect  of  the  (l+3x2)  * 

factor  is  to  increase  the  effective  wavelength  when  x  4  0.  Of  course, 
the  separability  evident  in  the  solutions  in  the  absence  of  this 
factor  disappears. 

One  might  expect  that  on  each  line  of  constaut  x  that  the 
radial  dependence  of  a  single  term  of  (3-25)  might  still  be  picked 

out.  This  is  the  case  but  the  value  of  m  is  no  longer  such  that  o< 

in 

on  that  particular  line,  is  a  minimum.  For  example,  at  x  ■  -.6  we 

expect  m  ■  15  (  real)  or  m  ■  16  (  ofjy  Imaginary)  but  find 

instead  m  ■  14(13  internal  nodes).  The  situation  is  best  demonstrated 
at  x  ■  0  where  if  the  concept  of  the  effective  value  of  X  on 

l+3x2 

a  line  of  constant  x  were  valid  we  would  expect  the  two  curves  of 
Figure  (5-16)  to  coincide.  In  actuality,  the  l+3x  curve  has  m  ■  13 
while  the  constant  X  curve  has  n  ■  12. 

Finally,  we  consider  the  case  X  ■  X  ,  where  X  is  a 

l+3x2 


•ora 


640  radial  points 


640  radial  points 


r  =  8,153  km 
20  angular  points 


r  =  16,963  km 
20  angular  points 


angular  points 


function  of  r.  X^  is  computer!  on  the  basis  of  the  model  ionosphere 
of  Prince  and  Bostick^^with  os  ■  1  cps. 

Figures  (5-21)  and  (5**22)  show  Xq  for  daytime  and  nighttime 
sunspot  minimum,  respectively. 

Figures  (5-23)  and  (5-24)  show  the  radial  and  angular 

functions  for  a  small  cavity,  r 2  ■  8,370  km.  Xq  is  plotted  next 

to  the  radial  function  and  it  is  clear,  for  this  particular  case, 

that  the  radial  function  is  insensitive  to  the  detailed  structure  of  X  . 

o 

Figures  (5-25)  and  (5-26)  show  the  radial  and  angular  func¬ 
tions  for  a  large  cavity,  -  19,880  km  for  the  nighttime  case. 

Figures  (5-27)  and  (5-28)  show  the  analogous  daytime  case.  Xq  is 
larger  in  the  daytime  and  hence  the  wavelength  is  smaller.  This 
produces  an  additional  node  in  the  radial  function  in  the  daytime  case. 
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x!09. 03  km 


5  angular  points 


6.  Conclusion 


This  report  has  described  the  successful  development  of  a 
computer  program  for  solving  the  poloidal  (unguided)  magnetohydrodynamic 
wave  equation.  The  solution  to  this  problem  is  very  closely  related 
to  the  determination  of  the  field  amplitudes  in  the  complex  asymmetric 
model  describing  the  real  ionosphere.  The  reason  for  this  is  that 
recent  developments  have  shown  that  neither  the  toroidal  (guided) 
mode  nor  the  highly  asymmetric  poloidal  modes  are  viable  „'n  a  dipole 
magnetic  field. 

We  have  investigated  the  credibility  of  this  computer  program 
by  comparing  numerical  results  with  analytic  results,  where  possible, 
and  with  the  results  of  a  different  computational  method. 

Thus  far,  there  has  not  appeared  to  be  any  intrinsic  difficulty 
in  incorporating  features  o_  the  real  ionosphere  into  the  model  although 
we  have  not  begun  to  explore  asymmetric  situations,  such  as  the  day- 
night  asymmetry. 
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APPENDIX 


TWO-POINT  BOUNDARY  PROBLEM 


The  two-point  boundary  problem  may  be  stated  as  follows: 

For  a  second  order  differential  equation  L  y(x)  =  S  where  y(x)  is 
the  unknown  function,  L  is  a  differential  operator  and  S  is  an 
inhomogeneous  term,  the  boundary  conditions  are  that  f  (y(x  ))  and 
f2(y<xb>>  are  specified  where  f^  and  f2  are  arbitrary  functions  or 
operators.  The  simplest  boundar3’  conditions  of  this  type,  which  are 
the  ones  we  are  actually  concerned  with,  are  y(xa>=0  and  yCx^^O. 
These  conditions  would  apply,  for  example,  to  a  stretched  string 
which  is  held  fixed  at  two  points,  x^  and  x^,  under  a  driving  force 
S.  If  the  driving  force  is  not  present,  then  such  a  problem  has  a 
non-trivial  solution  only  when  a  specified  parameter  has  certarn 
fixed  values,  i.e.,  eigenvalues.  When  there  is  a  driving  term  a 
solution  usually  exists  for  any  value  of  the  specified  parameter. 
However,  the  solution  is  usually  largest  when  the  parameter  equals 
one  of  the  eigenvalues  (i.e.  resonance  occurs).  Since  we  have  a 
driving  term,  the  problem  is  not  to  determine  the  eigenvalues  and 
eigenfunctions,  but  rather  to  solve  Ly=S  for  y  subject  to  the  two- 
point  boundary  conditions. 

In  order  to  solve  the  equation  numerically,  there  are  two 
standard  methods:  (1)  Total  solution  method;  (2)  Shooting  method. 

1.  Total  Solution  Method 

In  the  total  solution  method,  the  differential  equation  is 
differenced  at  a  set  of  grid  points  including  the  endpoints  so  that 
the  boundary  conditions  are  incorporated  directly  into  the  difference 
equations.  The  set  of  coupled  equations  leads  to  a  problem  of  solv¬ 
ing  a  set  of  linear  equations  involving  a  sparsely  populated  matrix 
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whose  order  is  equal  to  the  number  of  grid  points.  For  good  accuracy 
(small  truncation  error)  this  may  involve  using  a  large  number  of 
points.  Thj.s  method  was  used  in  the  case  where  the  partial  differ¬ 
ential  equation  (4-4)  was  differenced  in  two  dimensions.  It  has 
several  advantages  over  the  method  used  below,  but  unfortunately, 
as  of  this  writing  we  could  not  get  the  program  to  run  on  the  IBM 
7094  although  it  had  been  working  on  the  CDC  6400.  (There  seem  to 

be  some  limitations  in  the  former  machine  which  do  not  exist  in  the 
latter. ) 


2.  Shooting  Method 

Assuming  that  y(xfl)  is  given,  it  is  necessary  to  know 
y  (xfl)  in  order  to  start  the  numerical  integration.  In  the  shooting 
method  a  value  for  y'(xfi)  is  guessed  and  the  integration  is  performed. 
Then  the  value  y^),  if  this  is  given  as  the  boundary  condition  at 
xb,  would  be  compared  with  the  value  y(xb)  computed  on  the  basis  of 
the  assumed  value  of  y’(xa).  A  second  integration  would  be  per¬ 
formed  with  a  different  value  of  y’(xa).  By  suitably  choosing 
y'(xfl)  the  solution  can  be  made  to  assume  the  value  y^)  specified 
by  the  boundary  conditions. 


This  method  has  several  disadvantages:  (1)  It  is  not  an 
exact  method.  Several  integrations  must  be  performed  before  a 
satisfactory  value  for  y'(xfl)  is  found,  that  is,  before  a  solution 
is  found  whose  value  y(xfe)  is  within  a  prescribed  tolerance.  Also, 
and  probably  more  significant,  if  there  is  any  tendency  toward 
instability,  it  will  be  exhibited  when  the  equation  is  being 


integrated  for  incorrect  values  of  yf(x  ). 

a 

In  spite  of  these  difficulties  it  is  the  second  method, 

the  shooting  method,  which  has  been  used  successfully  i:  the  present 
case. 
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For  the  case  when  Eq.  (4-4)  is  tc  oe  solved,  this  is 
equivalent  to  solving  2N-2  first  order  differential  equations 


7; .  l 


IN -2  ) 


) 


(A— 1) 


Here  the  y  with  odd  subscripts  are  the  amplitudes  and  y  with  even 
subscripts  are  the  derivatives.  The  boundary  conditions  are  that 
y^Cx^-O  and  yi(x^)=0  for  i  odd.  The  assumed  values  for 

are  denoted  by  n  .  Since  the  set  (n.)  is  chosen  anew  for  each 

1  1  0) 
integration,  this  will  be  denoted  by  a  superscript,  i.e.  q  * 

where  1  <  j  <  N.  The  actual  values  achieved  for  y0^  .^(x^)  w^en 

are  used  as  inputs  are  denoted  by  .  The  assumption  is  made  that 


£ 


(A-2) 


where  a^  and  b^  are  constants.  That  is,  it  is  assumed  that  there 

is  a  linear  relationship  between  the  set  {ru  }  and  the  set  {? }. 

(1)  K  1 

Then  we  have,  using  the  quantities  ,  found  by  integrating  the 
differential  equations 


N-l 

X 


(A-3) 
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Now  for  each  i  the  N  equations  (i.e.,  1  <  j  <  N)  are  a 

linear  set  of  equations  for  the  unknowns  a^  and  b^.  Solving  these 

we  find  the  a.,  and  b.  so  that  we  have 
ik  i 


f 


(J-t 


(A-4) 


If  we  now  set  £ 


0  as  demanded  by  the  condition  =  0 


(i  odd),  then  the  'Jj  are  determined  by  solving 


/V-/ 


(A-5) 


This  gives  a  new  set  of  jf-  from  which  to  start  the  integration  and 
this  leads  to  a  new  set  of  . 

The  next  step  in  the  method  is  to  test  the  <5^  to  see  if 
they  are  within  acceptable  bounds.  If  so,  the  solution  has  been 
found.  If  not,  a  new  calculation  is  performed  for  a  and  b 

-  ^  _ .  _  IK  i_ 

using  the  sets  of  <**«/  j  fi  J  for  2  <  j  <  N+l.  From  this  a  new 

set  of  TK  are  predicted,  integration  of  the  equations  performed  and 
the  new  are  tested.  This  process,  in  theory,  continues  until 

an  acceptable  solution  is  reached.  In  practice,  it  has  been  found 

»  -6 

impossible  to  make  the  individual  less  than  about  10  of  the 
derivatives  at  toe  end  point.  This  is  related  to  the  fact  that 
the  function  and  its  derivative  are  related  quantities  and  that 
there  are  about  8  significant  figures  held  on  the  IBM  machine.  When 
calculations  were  performed  on  a  CDC  machine  which  held  more  signifi¬ 
cant  figures  for  each  number,  the  convergence  was  improved. 


-59- 


One  interesting  facet  of  this  technique  is  related  to  a 
view  in  which  the  (tl^)  are  considered  imputs  to  a  system  (of 
differential  equations)  for  which  (  ?  )  are  the  outputs.  In  order 
to  solve  (A-3)  for  aik  we  must  deal  with  the  matrix  of  imputs 

which  we  are  free  to  choose. 

Since  the  matrix  is  not  to  be  ill-conditioned  or  singular 

we  require  that  the  numbers  not  be  chosen  in  ,my  simply-related 

manner  because  in  that  case  the  matrix  would  be  singular  or  close 

to  it.  In  particular,  we  have  chosen  the  elements  of  tif^  to  be 

k 

random,  so  that  the  method  used  here  to  sclve  the  differential 
equation-two  point  boundary  problem-may  be  viewed  as  equivalent 
to  exploring  a  system  by  means  of  observing  its  response  to  random 
imputs. 


B.  EXPERIMENTAL  STUDY  OF  ION- ACOUSTIC  WAVE  PROPAGATION 


1.  Introduc t  ion 


One  of  the  problems  in  the  magnetic  detection  of  a  high 
altitude  nuclear  detonation  is  in  predicting  the  propagation  of  the 
wave  signal  through  the  lower  ionosphere  where  the  plasma  and  neutral 
density  vary  rapidly  with  altitude.  The  problem  is  particularly 
difficult  for  an  analytical  approach  because  the  plasma  model  varies  in 
this  region  from  a  collisionless,  fully  ionized ' plasma  to  that  of  a 
collision-dominated,  weakly  ionized  gas.  The  wave  signal  which 
propagates  through  this  region  is  attenuated  by  collisions  with  the 
neutrals.  Such  collisional  effect  couples  the  plasma  wave  motion 
into  the  neutrals,  and  the  coupling  is  strong  when  the  collision 
frequency  becomes  comparable  with  the  wave  frequency.  Existing 
theoretical  studies  on  the  problem  are  based  on  various  simplifying 
models.  The  validity  of  these  models  can  only  be  assessed  by 
laboratory  experiments  under  controlled  conditions. 

The  laboratory  technique  for  studying  the  phenomena  of 
interest  has  been  discussed  in  a  previous  report^^ .  Efforts  during 
the  current  reporting  period  have  been  directed  toward  construction 
of  the  experimental  apparatus,  instrumentation,  and  experimental 
planning.  This  work  is  reported  in  the  following  sections. 
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2.  Construction  of  the  Experimental  Apparatus 


It  was  determined  that  the  experiment  can  best  be  carried 
out  in  an  alkaline  plasma  generator  operating  in  the  vapor  pressure 
mode.  In  order  to  generate  such  a  plasma,  several  pieces  of  equipment 
are  needed.  These  include  a  vacuum  manifold  and  associated  pumping 
network,  a  temperature  regulated  cesium  reservoir  for  maintaining 
the  cesium  vapor  pressure,  an  oven  to  maintain  the  entire  assembly 
at  an  elevated  temperature,  and,  finally,  an  electron  bombardment 
heating  system  for  heating  the  ionizing  hot  plate.  The  design  and 
construction  of  the  equipment  are  described  in  the  following  sections. 

A.  Vacuum  Manifold 

The  vacuum  manifold  serves  to  confine  the  alkaline  metal 
vapor  in  a  controlled  environment.  The  vacuum  system  is  shown 
schematically  in  Fig.  1.  The  system  consists  of  a  cylindrical 
vacuum  manifold  20  cm  in  diameter  by  100  cm  in  length,  enclosed 
by  blank  flanges  on  each  end.  Metal  gaskets  are  used  throughout  to 
insure  high  temperature  operation  of  the  system.  A  pumping  duct  is 
connected  to  the  system  through  one  of  the  flanges.  The  system  can 
be  isolated  from  the  pumps  during  the  experiment  by  the  bakeable 
Granville-Phillips  type  267  Auroseal  valve.  This  valve  was  chosen 
because  of  its  ability  to  operate  at  elevated  temperatures,  both  in 
the  open  and  closed  positions.  A  cold  trap  is  inserted  in  the  pump¬ 
ing  line  to  condense  out  the  room  temperature  cesium  vapor, which  can 
attack  the  pumps  and  also  present  a  health  hazard.  When  the  system 
is  isolated  from  the  pumping  network  after  initial  pump  down,  cesium 
vapor  is  introduced  into  the  system  by  heating  up  the  reservoir, which 
is  connected  to  the  system  via  a  Tee  in  the  pumping  duct.  The  tem¬ 
perature  of  the  reservoir  is  controlled  by  a  heating  element  imbedded 
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in  an  insulating  jacket  surrounding  it.  The  flanges  are  equipped 
with  vacuum  feedthru  elements  for  power  and  signal  connections  into 
the  plasma.  A  view  of  the  vacuum  system  assembly  is  shown  in  the 
photograph  in  Fig.  2. 

B.  Oven  Assembly 

In  order  to  avoid  the  condensation  of  the  cesium  vapor, 
it  is  necessary  to  maintain  the  entire  system  at  an  elevated  tempera¬ 
ture.  It  was  originally  intended  that  the  entire  system  be  maintained 
in  an  oven  heated  by  the  hot  plate.  However,  later  heat  transfer 
calculations  have  shown  that,  because  of  the  large  heat  flux  escaping 
from  the  source,  it  would  be  desirable  to  cool  the  regions  directly 
heated  by  the  source  to  avoid  excess  heating  and  separately  heat  the 
pumping  ducts  >md  the  bakeable  valve  by  external  heating  elements. 
Thus,  the  fiv  .1  design  consists  of  thermally  insulated  pumping  ducts 
and  valve  with  imbedded  heating  elements  similar  to  the  cesium 
reservoir,  and  the  entire  system  is  installed  on  a  rack  to  which 
cooling  fans  can  be  directed  to  provide  cooling  desired. 

C.  Electron  Bombardment  Heating  System 

The  cesium  plasma  is  to  be  generated  by  surface  ionization 
of  the  incotfiing  neutrals.  In  order  to  maintain  charge  neutrality, 
electrons  must  be  evaporated  off  as  rapidly  as  ionization  takes  place. 
The  electrons  are  to  be  generated  by  thermionic  emission.  The 
temperature  of  the  hot  plate  required  to  maintain  equal  electron 
emission  and  ionization  is  shown  in  Fig.  3  as  a  function  of  neutral 
temperature  for  cesium  on  tungsten. 

To  assure  that  the  plasma  sees  only  a  well  defined  equi- 
potential  surface,  it  is  desirable  that  the  hot  plate  be  indirectly 
heated.  Electron  bombardment  heating  is  ideally  suited  for  that 
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purpose.  In  order  to  generate  a  plasma  of  sufficient  volume,  it 
was  determined  that  a  surface  area  of  5  cm  would  be  required.  A 
zig-zag  filament  covering  a  circular  area  of  10  cm  diameter  was 
formed  with  tungsten  powder.  The  filament  is  to  be  heated  by  a  10- 
volt,  200-amp, current  regulated  power  supply.  The  filament  is 
biased  at  -5000  volts,  and,  thus,  the  thermionically  emitted  electrons 
are  accelerated  to  the  tungsten  hot  plate  where  the  kinetic  energy 
is  dissipated  in  heating  up  the  plate.  In  order  to  attain  uniform 
temperature  on  the  surface  of  the  hot  plate,  a  relatively  heavy 
plate  of  1  cm  thickness  was  used. 

To  heat  the  plate  to  a  high  temperature,  a  rather  sub¬ 
stantial  power  supply  with  good  regulation  characteristics  is 
required  to  maintain  operational  stability.  The  high  current 
required  to  heat  the  filament  and  the  high  voltage  required  to 
accelerate  the  electrons  are  incorporated  into  a  single-unit  design 
which  produces  a  maximum  current  stabilized  supply  of  200  amps  at 
10  volts,  or  100  amps  at  20  volts  to  the  filament  biased  at  -5000 
volts.  The  high  voltage  supply  is  voltage  regulated  and  continu¬ 
ously  variable  from  500  to  5000  volts  with  a  maximum  output  of 
1.5  amps.  The  construction  of  the  power  supply  has  been  completed 
and  the  unit  is  now  operational. 

The  three  items  discussed  above  represen.  the  major 
equipment  required  for  the  experiment,  and  construction  of  all  of 
these  items  has  been  completed. 
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3.  Instrumentation 


Current  work  In  the  progress  of  the  experiment  has  been 
concentrated  in  instrumentation  of  the  apparatus.  This  instrumenta¬ 
tion  includes  the  installation  of  the  hot  plate  assembly,  diagnostic 
sensors,  and  radiation  shields.  Because  the  system  is  to  be  operated 
at  an  elevated  temperature,  metal  gaskets  are  required  for  all 
demountable  seals.  Therefore,  it  was  decided  at  an  early  stage,  for 
economic  reasons  as  well  as  for  simplicity,  that  all  instrumentation 
is  to  be  Installed  inside  the  vacuum  manifold  with  all  leads  brought 
out  through  feedthru  elements  installed  in  the  end  flanges.  To 
support  the  /arlous  electrodes ,  two  rods  cantilevered  at  one  of  the 
end  flanges  run  the  length  of  the  vacuum  manifold.  The  hot  plate 
assembly,  electrostatic  probes  and  grids,  thermocouples,  and  radiation 
shields  are  all  hung  from  these  two  rods.  The  detail  construction 
of  the  various  parts  is  described  in  the  following. 

A.  Hot  Plate  Assembly 

The  hot  plate  consists  of  a  1  cm  thick  tungsten  circular 
plate,  10  cm  in  diameter,  mounted  in  two  boron  nitride  disks  as 
shown  in  Fig.  4.  The  boron  nitride  disks  serve  as  electrical  insula¬ 
tion  so  that  the  hot  plate  bias  may  be  varied  with  respect  to 
ground.  The  disks  are  held  together  by  two  brackets  made  from 
tantalum  sheet  metal  and  attached  to  the  supporting  rod.  The 
tungsten  filament  is  attached  directly  onto  the  high  current,  high 
voltage  feedthru  elements,  and  the  entire  assembly  is  surrounded  by 
radiation  shields  to  cut  heat  loss.  This  design,  differs  from  the 
one  proposed  in  the  last  semi-annual  report  in  greater  simplicity 
and  ease  of  installation. 
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AT  ION  SHIELDS 


B.  Diagnostics 


The  basic  instrumentation  for  the  experiment  consists  of 
electrostatic  probes  for  measuring  spaiial  distribution  of  the  plasma, 
grids  for  excitation  and  detection  of  waves,  and  thermocouples  to 
monitor  the  temperature  in  the  reservoir  and  other  regions  of  interest. 
The  probes  and  grids  are  fabricated  from  tungsten  wires  insulated  by 
alumina  tubing  and  mounted  on  a  specially  designed  probe  holder  which 
will  self-center  the  probe  at  any  radial  and  axial  positions  desired. 
The  signals  are  transmitted  by  Teflon  insulated  wires  which  can 
withstand  temperatures  of  500°K. 

The  temperature  of  the  reservoir  is  monitored  by  a  thermo¬ 
couple  attached  to  the  outside  metal  wall.  Since  the  reservoir  is 
enclosed  in  heavy  insulation,  we  anticipate  the  temperature  drop 
across  the  metal  wall  will  be  insignificant,  and  thus  the  external 
temperature  measurement  will  give  a  satisfactory  estimate  of  the 
reservoir  temperature.  The  neutral  b.  Aground  pressure  in  the  system 
will  be  ascertained  from  the  measured  temperature  of  the  reservoir 
and  the  vapor  pressure  curve  of  cesium. 

To  avoid  condensation  of  the  vapor  at  any  point  ir.  the 
system,  it  is  important  to  assure  that  the  temperature  throughout  the 
system  is  kept  above  that  of  the  reservoir.  Since  the  only  require¬ 
ment  for  temperature  in  this  case  is  that  it  should  be  above  a  certain 
value,  we  have  painted  the  system  at  various  control  points  with 
Thermolaq  paint  rather  than  actually  measuring  temperatures  at  these 
points.  This  procedure  considerably  simplified  the  monitoring 
requirements. 


C.  Radiation  Shields 


In  order  to  maintain  the  hot  plate  at  elevated  temperatur 
with  minimum  power  input,  and  also  to  protect  the  wire  insulations 
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from  over  heating,  it  is  desirable  to  minimize  the  heat  transfer 
from  the  central  regions.  Calculations  have  shown  that  radiation 
and  free  convection  are  able  to  remove  about  500  watts  from  the 
vacuum  walls  at  500°K.  Forced  convection  can  increase  this  rate, 
possibly  by  a  factor  of  two.  Therefore,  radiant  heat  transfer  from 
the  hot  plate  must  not  exceed  this  value.  This  can  be  accomplished 
by  the  installation  of  10  layers  of  radiation  shields.  Present  plan 
calls  for  the  installation  of  15  layers  of  1  mil  thick  tantalum  foil 
surrounding  a  useful  volume  of  about  3  litres.  This  shield  is 
presently  being  installed. 

The  final  experimental  work  can  be  carried  out  upon  the 
completion  of  these  various  installations  in  the  near  future. 


4.  Experimental  Program 


The  basic  objective  of  the  experiment  is  to  study  the 
propagation  of  ion  acoustic  waves  in  a  plasma  medium  where  the 
parameters  may  be  varied  under  controlled  conditions.  The  basic 
experimental  arrangement  is  shown  schematically  in  Fig.  5.  Longi¬ 
tudinal  waves  are  excited  in  the  plasma  by  applying  a  tone  burst 
of  a  chosen  frequency  oscillation  and  detected  by  sensors  along  the 
plasma  column.  A  tone  burst  signal  is  preferred  over  a  CW  signal 
so  that  directly  coupled  signals  can  readily  be  distinguished  from 
the  ion  wave  signals. 

The  controlled  experimental  variables  are  the  plasma 
density  n,  the  neutral  density  nQ,  the  wave  period  t,  the  scale 
length  and  grid  separation  d^.  The  experimentally  measurable 
parameters  of  interest  are  the  time  delay  ^  and  of  the  signal 
arriving  at  detectors  1  and  2  and  the  amplitude  of  the  received 
signals  and  A^.  From  these  data  one  can  construct  the  dispersion 
relation  for  the  system  which,  in  principle,  defines  the  system 
completely. 

Comparison  of  the  experimentally  determined  wave  propagation 
characteristics  with  theoretical  predictions  can  be  expected  to  yield 
valuable  insight  into  theoretical  modelling  of  waves  in  inhomogeneous 
media.  The  phenomenon  of  particular  interest  is  the  damping  of  the 
waves.  It  is  expected  that  we  will  be  able  to  observe  both  collisional 
and  collisionless  damping  which  are  manifestations  of  various  coupling 
mechanisms.  The  optimal  detection  method  for  signals  propagating 
through  such  a  weakly  ionized  medium  can  be  determined  from  the 
experimental  conclusions. 
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